Intermittent quakes and record dynamics in the thermoremanent magnetization of a 

spin-glass 
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A method for analyzing the intermittent behavior of linear response data in aging systems is 
presented and applied to the spin-glass thermoremanent magnetization (TRM) data of Rodriguez 
et al. (Phys. Rev. Lett. 91, 037203, 2003). The probability density function (PDF) of the 
magnetic fluctuations has an asymmetric exponential tail, showing that the demagnetization process 
occurs through intermittent spin rearrangements or quakes which significantly differ from reversible 
fluctuations having a Gaussian distribution with zero average. The intensity of quakes is determined 
by the TRM decay rate, which in turn depends on t, the time since the initial quench and on t w , 
the time at which the magnetic field is cut. For a broad range of temperatures, these dependences 
are extracted numerically from the data and described analytically using the assumption that the 
system's linear response is fully subordinated to the occurrence of the quakes which spasmodically 
release the imbalances created by the initial quench. 

PACS numbers: 65.60.+a, 05.40.-a, 75.10.Nr 
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INTRODUCTION 

Many glassy materials unable to re-equilibrate after 
a temperature quench undergo aging, a dynamical pro- 
cess whose key properties are largely system independent. 
In spin-glasses, a class of disordered magnetic materials, 
one thoroughly investigated aging quantity is the ther- 
moremanent magnetization (TRM) [l|. In a TRM ex- 
periment, a sample undergoes a rapid thermal quench at 
t = 0, with a small magnetic field present. Subsequently, 
the temperature remains constant throughout the exper- 
iment, while the magnetic field is turned off in a single 
step at t = t w . The decay of the TRM for t > t w de- 
pends mainly on the scaling variable (t — t w )/tl^, where 
fi is a parameter close to unity [l|- As recently recog- 
nized 2,3], the value of fi depends on the quenching rate, 
and approaches a so-called full aging limit, /i = 1, as the 
rate increases toward infinity, i.e. in the limit of an in- 
stantaneous quench. These phenomena together reveal a 
persistent memory of the initial quench and a strong sen- 
sitivity to the rate of cooling. Memory behavior includes 
a number of other fascinating aspects 0, [B| , and could be 
rooted in a multi-scale, hierarchical nature of the energy 
landscape of amorphous systems [1, 0, H, @, EU • 

Memory issues can be further elucidated by analy zing 
the fluctuations statistics in meso-scaled systems [lj, QJl, 



13l . Il4j . The basic observation is that the Probability 



Density Function (PDF) of an aging quantity typically 
comprises a Gaussian part and an exponential tail. The 
Gaussian covers reversible fluctuations with zero average, 
and the tail describes intermittent events, which carry 
the drift of the aging process, e.g., in the present case, 
they carry the net change in magnetization. The inter- 



mittent events are prominent during the non-equilibrium 
aging regime t » t w where the Fluctuation Dissipation 
theorem is clearly violated 12, [3, EH- Furthermore, 
as shown in the sequel, they are equally prominent in a 
short time interval immediately following field removal. 
In general, the statistical weight of intermittent fluctu- 
ations relative to the weight of Gaussian fluctuations is 
determined by the rate ttrm of magnetization change, 
the latter depending on both t and t w as detailed below. 

Heat transfer statistics directly probes thermally acti- 
vated dynamics, and is easily collected in numerical sim- 
ulation. In two different agingsystems, the Edwards- 
Anderson spin-glass model [a, LL5| and a model with p- 
spin interactions (lij . Gaussian and intermittent events 
are clearly identifiable. The idea that the intermittent 
events, or quakes are triggered by extremal, e.g. record 
sized, reversible fluctuations 17|, [l8( explains the age de- 
pendence of the PDF of energy fluctuations. The same 
'record dynamics' approach was applied to the config- 
uration autocorrelation function [191 ] and to dynamical 
properties of other complex systems [20|, |2l| . 

Arguably, quakes influence aging quantities other than 
the energy, including e.g. the all important linear re- 
sponse functions. Unlike the energy, response functions 
involve the turning of an external perturbation on or off, 
and hence require additional theoretical attention. In 
equilibrium situations, it is well known that a small per- 
turbing field simply probes the spontaneous fluctuations 
of its conjugate variable. For out-of-equilibrium intermit- 
tency, an idea similar in spirit can be formulated as a sub- 
ordination hypothesis: significant configuration changes, 
whether field induced, as in TRM, or spontaneous, as in 
the de-correlation of the magnetization fluctuations of an 



2 



unperturbed system, only occur in conjunction with the 
quakes. These, in turn, release the strain created by the 
initial quench in a temporal sequence which is, to linear 
order, unaffected by the perturbation. This subordina- 
tion hypothesis was recently applied to the autocorrela- 
tion function in the Edwards- Anderson spin-glass in zero 
field [ijj and is presently applied to the experimental 
TRM decay. At a price — short time effects of equilib- 
rium like-fluctuations are neglected — the approach leads 
to a considerable mathematical simplification, as the sta- 
tistical properties of e.g. the response flow from those of 
the quakes. 

In the sequel, we first describe the method of data 
analysis. We then show the intermittency of the decay 
and fit the age and time dependence of the TRM decay 
rate to a simple formula. The theoretical ideas behind 
this formula are spelled out in the Theory Section, and 
the whole paper is rounded off with a brief concluding 
Section. 

Last but not least, a notational issue: the variable V 
denotes here the time elapsed since the initial quench, 
or 'age', and 't w ' denotes the specific age at which the 
external field is switched off. A slightly different con- 
vention, widespread in the experimental literature, see 
e.g. Refs. [H, 0, [53|, uses 'f for the time elapsed after 
field removal, a quantity presently denoted by 'i &s\ with 



tobs =t — t w . Furthermore, in refs. [151. 124j. t w is used to 
denote the system age, our present t. Our choice em- 
phasizes the initial quench as the common origo for all 
time variables, and slightly simplifies the notation, e.g. 
our scaling functions have a t/t w rather than a 1 + t/t w 
argument. 



METHOD OF DATA ANALYSIS 

In this Section, the TRM data of ref. || are ana- 
lyzed with focus on intermittency. In the experiment, 
a Cuo.94Mno.06 sample, with critical temperature T g = 
31. 5K, is rapidly cooled, with an effective cooling time 
of 19s. The sample is then aged isothermally at tem- 
perature T in a weak magnetic field (H = 20G), with 
the field removed at t = t w . The TRM signal M(t) 
is measured for t > t w and for different values of t w 
and T: t w = 50, 100, 300, 630, 3600, 6310 and 10000s and 
T = 0.4,0.6,0.83,0.9 and 0.95 T g . The observation time 
tobs = t — t w ranges from appr. 6000s to appr. 30000s 
and the ratio t/t w correspondingly varies from 1 to 100. 
Measurements are recorded every 1.045s for t w < 6310s 
and every 2.04s for t w > 6310s. All data are given as 
dimensionless ratios of the TRM to the field cooled mag- 
netization Mpc- 

We consider TRM changes SM over a small time in- 
terval St, i.e. from the data we calculate a time series of 
magnetization differences 8M(i) = M(t$ + St) — M(tj), 
with ti = t,j_i 4- St, and with St << t chosen as a small 



multiplum of the measurement repeat time. The analysis 
has a twofold aim: estimating the PDF of the magnetiza- 
tion fluctuations in different situations in order to show 
the presence of intermittency, and estimating the time 
and age dependence of the rate ttrm (henceforth simply 
'rate') of magnetization decay, in order to compare with 
theory. 

The PDF are straightforwardly estimated by binning 
the 5M(i) values sampled over suitable time intervals. 
An estimate of ttrm is not easily obtained: in order 
to uncover the drift part of the dynamics, the Gaussian 
fluctuations, which are far more frequent than intermit- 
tent events over small time intervals, must be averaged 
out. Mainly, the averaging is done over suitably chosen 
subintervals of the observation interval. Additionally, we 
perform for T — 0.83T g an average over an ensemble 
consisting of data with the same t/t w value. From this 
ensemble we also construct the PDF of the magnetic fluc- 
tuations over a short time interval near t w . 

For time averaging, we consider a set of intervals It = 
[(8/9)£, (10/9)f] with midpoints t equidistantly spaced on 
a logarithmic scale. The interval widths, 0.22i, are chosen 
as a compromise between the conflicting requirements of 
small statistical error and good temporal resolution. Us- 
ing the values available in each interval I t , the average 
USAdit, t w , St) and the variance cr^ M (i, t w , St) of the mag- 
netization change SM are estimated at time t using stan- 
dard formulas. By varying St, we ascertain that fj,$M is 
proportional to St. The proportionality constant, which 
is identical to rTRM{t,t w ), is estimated by linear regres- 
sion. The t and t w dependence of ttrm thus obtained is 
fitted to the expression predicted by record dynamics. As 
a final step, the dependence is integrated with respect to 
t, leading to an analytical expression for the TRM decay, 
which is compared with the original data. The outcome 
of this whole procedure is displayed in the first five panels 
of Fig. El 

For ensemble averaging, we note that the t/t w depen- 
dence of the TRM (full aging), which is established in the 
literature Q and confirmed by the present analysis, im- 
plies that SM values with the same t/t w , i.e. t/t w = C, 
are physically equivalent. Hence, they can meaningfully 
be collected from data streams taken at different t w into 
ensembles labeled by the value of C. This procedure is 
carried out for the T = 0.83T g data, collecting, in each 
case, SM values within the interval t £ [Ct w ,1.25Ct w ] 
and systematically varying the value of C. From the 
C = 1.13 ensemble, we construct the (unnormalized) 
PDF of SM data. We also estimate, for a number of dif- 
ferent C, the rate rxRM{t,t w ) as the average of SM/St 
of the corresponding ensemble. The results are shown in 
the main panel and the insert of Fig. [31 respectively. 
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(a)Intermittent and Gaussian PDF (b)Average and 



FIG. 1: (Color on line) (a): The circles show the PDF of the magnetization change SM for isothermal aging at T — 0.83T 9 , 
with the field cut at t w = 100s. The values of SM are taken over short intervals St = 1.045s within the observation interval 
[1000, 4000]. Diamonds show the PDF of SM from the same data stream, but now taken over larger intervals St = 3 x 1.045s. 
This leads to a much stronger intermittent component on the left wing of the PDF. The full lines are least square fits of the 
positive values of SM to a Gaussian with zero average, (b): From the same data, the average and the variance (the latter scaled 
as shown) of SM's are plotted versus St. The dotted line is obtained by a least square error fit. 



RESULTS 

Figure [1] illustrates how TRM fluctuations SM occur- 
ring over a small interval St have both a Gaussian com- 
ponent and an intermittent tail. The tail becomes more 
dominant as the average change lism grows numerically 
larger. This differs from 'normal' transport, where in- 
creasing fisM would only shift the center of the Gaussian 
distribution away from zero. 

Panel (a) of Fig. Q] compares two (unnormalized) PDFs 
of SM, which are both obtained from data with T — 
0.83T g and t w = 100s. The SM(i) values are collected 
over the same observation interval / = [1000s, 4000s]. 
The PDFs are shown on a logarithmic scale, where a 
Gaussian has a parabolic shape. The nearly Gaussian 
PDF (circles) is for St = 1.045s and the other PDF (di- 
amonds) is for St = 5 x 1.045s. Increasing St increases 
the average magnetization change fiSM — St ttrm and 
hence increases the amount of intermittency. Using pos- 
itive SM values, the central part of the PDF is fitted 
(full lines) to a Gaussian with zero average. Note how 
the Gaussian shape of the reversible fluctuations remains 
visible for positive SM values, in spite of the strong in- 
termittent left wing present for the larger St (diamonds). 

Panel (b) of Fig. Q] shows the average fisM (squares) 
and the variance a$ M (diamonds), plotted versus St. As 
fiSM and St are expectedly proportional, the rate of mag- 
netization decay, (averaged over the observation interval 
I) can be estimated as the slope ttrm = H8Ai/St. The 
dotted line through the data illustrates the quality of the 
linear fit. For completeness, we also show the variance 



Cg M . As one would expect in a large sample, agM is 
much smaller than the average. 

Panels 1-5 of Fig[2]illustrate the t and t w dependence of 
the rate of magnetization ttrm for the five temperatures 
indicated. Error bars are obtained by standard methods 
and data points with (la) uncertainty larger than 10% 
are discarded. As it turns out, the rate approaches the 
functional form ttrm k a /t for large values of t/t w , i.e. 
already for t/t w > 10 this is the main contribution to 
Ttrm- Thus, t Ttrm — a plotted versus t/t w describes 
the deviation of the rate from its asymptotic behavior. 
The full line depicts the function 



With parameters optimized by a least square error fit, 
this function offers a good analytical description of 
t ttrm —o,- The form of the parameterization is justified 
theoretically in the next Section. Here we note that as 
A 2 < —3, the second term of Eq. [T]only contributes for 
t » t w , and that the sole parameter of importance in the 
limit t » t w remains a, which, importantly, remains 
nearly constant through the temperature range. 

Let us finally consider the change in magnetization 
AM over an observation interval I — [ti,t] starting at 
a time U larger or equal to t w , but otherwise arbitrary. 
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Plainly, 

AM(ti,t,t w ) 



Hi 

h 
Ai 

A 2 



rTRM{t' \t w )dt' = a\og(t/t % ) 

Ai / , \ Ai" 



i"tt 



(2) 



For a generic observation interval, AM is a function of 
three variables, as indicated. Customarily, one chooses 
ti = t w , and indeed, for ti — t w — 100s, the above for- 
mula yields the analytical approximation AM(t,t w ) to 
the TRM decay plotted (blue circles) for selected values 
of t, together with the measured data (red line) in the 
inserts of Fig. [2l Conforming to standard usage, the ab- 
scissa t/t w — 1 is the ratio of the observation time to the 
field removal age. Note that the t w dependent value of 
the magnetization at t = t w , formally an integration con- 
stant, is not provided by Eq. [5J This value is determined 
by shifting AM(t,t w ) vertically, until they best overlap 
with the data is obtained. Finally, the asymptotic de- 
cay of the magnetization only acquires a t w dependence 
if ti = t w is chosen. 

The last panel of Fig. [2] shows the temperature depen- 
dence of the parameters a, Ai and A2. Remarkably for 
thermally activated dynamics, the value a = —0.01 works 
for all temperatures except for T — 0.95T g , i.e. very close 
to the critical temperature. This behavior has a simple 
interpretation: for t/t w sufficiently large, each quake con- 
tributes, on average, the same amount to the TRM decay, 
and a is proportional to the logarithmic rate of quakes, 
a, which is T independent according to the theory. Fi- 
nally, the temperature dependence of the two exponents 
Ai and A2 and of the corresponding pre-factors b\ and 62 
(not shown), is smooth and rather weak, except near T g . 

A property not previously noticed is the strong inter- 
mittency for t rs t w . This is seen in the main plot of 
Fig.O where the (unnormalized) PDF (dots) of t5M/6t- 
a is collected from intervals [t, 1.25t] with t/t w = 1.13, 
using all available data streams with T = 0.83T g . The 
denominator St is the repeat time of the measurement. 
The PDF features a Gaussian component and a strong 
intermittent component, as highlighted by fits (full lines) 
to a Gaussian with zero average (obtained from the posi- 
tive data values) , and an exponential (obtained from the 
negative data values). As a consistency check, we esti- 
mate the quantity t Ttrm — a as a function oit/t w . The 
insert shows this quantity (dots) with ttrm estimated for 
each t/t w as the average of 8M/5t over the correspond- 
ing ensemble. The line — taken from the third panel of 
Fig. [2] — is a fit obtained via time averages. The agree- 
ment shown between time and ensemble averages (except 
at very small values of t/t w ) confirms the validity of the 
time averaging procedure used to evaluate ttrm- 



Summarizing, the fast magnetization decay occurring 
immediately after field removal is intermittent, a further 
indication that all magnetization decay is intermittent 
and controlled by the magnitude of ttrm- Secondly, 
time and ensemble averages give similar estimates of the 
magnetization rate. 



THEORY 

As it transpires from our results, the main theoreti- 
cal focus will be on the rate of demagnetization ttrm, 
which appears intimately related to the intermittent 
events. 



In a record dynamics scenario [171 lla ] , intermit- 
tent events, or quakes, reflect significant configurational 
changes which (i) lead from one metastable configuration 
to another, (ii) are irreversible and (in) are triggered by 
extremal fluctuations. For thermal activated dynamics, 
these would be thermal fluctuations over record-sized en- 
ergy barriers. We presently assume that quakes are sta- 
tistically independent, and that measurable effects are all 
subordinated [l9| to their occurrence. The time and age 
dependence of any aging quantity is then co-determined 
by how many quakes occur between t w and t and by the 
magnitude and nature of the physical changes that each 
quake entails. The number of quakes nj in the observa- 
tion interval, (t w ,t), is a Poisson distributed stochastic 



variable [171 1 1 Si ] with average 



(ni(t w ,t)) = a log(t/t w ). (3) 

The parameter a depends linearly on the system size [l^ | , 
as it arises from an extensive number of independent 
contributions from different thermalized domains. Fur- 



thermore, it is temperature independent [15[, due to the 
implied self-similarity of the energy landscape of each 
thermalized domain _[16| • While the above properties are 
supposedly generic [17| , the physical effects of the quakes 
could be system and even variable dependent. However, 
even without a specific knowledge of these effects, a sub- 
ordination hypothesis significantly restricts the possible 
age dependencies [lj|. of physical quantities of interest. 

Subordination means that any time dependence is me- 
diated by ni {t w ! t) , whence the TRM magnetization be- 
comes an (unknown) stochastic process with nj acting 
as an effective 'time' variable. In principle, desired prop- 
erties of the TRM can be found for fixed nj, and the 
underlying time dependencies can be re-introduced, by 
averaging nj according to the Poisson distribution spec- 
ified by Eq. El 

As quakes are seldom events, we can treat their phys- 
ical effects, e.g. magnetization or energy changes, as 
statistically independent. If the thermorcmancnt mag- 
netization is treated as a Markov chain with ni in the 
role of (discrete) time, its average Mt_rm(«/) admits an 
eigenvalue expansion [25( . According to Eq. [31 the range 
of nj will be modest for achievable time arguments and 
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FIG. 2: (Color on line) Panels 1-5: the main plots show the estimated values of t Ttrm — a versus t/t w , for t w = 50 (right 
pointing triangles, unavailable for T = 0.6T 9 and T = 0.9T 9 ), 100 (circles), 300 (squares), 630 (diamonds), 1000 (pentagrams), 
3600 (hexagrams), 6310 (asterisks) and (10000) (left pointing triangles). The lines are given by Eq. [T] The inserts compare 
the TRM decay measured at t w — 100s (red line), with the theoretical estimates (blue circles) obtained by integrating the 
fitted decay rate, see Eq. [2] The abscissa t/t w — 1 is the ratio of the observation time to the field removal age. In panel 6, the 
parameter a and the two exponents Ai and A2 are plotted vs. temperature. The lines are guides to the eye. 
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g 




t (5M/8t)-a 

FIG. 3: (Color on line) The main plot shows the unnormal- 
ized PDF (circles) of the quantity SM/St — a. The values 
of SM are collected from data streams taken with different 
t w and with T = 0.83T g . The denominator St is the repeat 
time in each measurement, All data are taken within inter- 
vals [t, 1.25£], where the value of t changes with t w , keeping 
a constant ratio t/t w — 1.13. The parabola and the straight 
line are fits to a Gaussian and an exponential, respectively. 
The insert shows the t/t w dependence of t TTRMt — a (dots), 
where ttrm is estimated by averaging SM/St over the en- 
semble of data points available at each t/tw The full line, 
which is lifted from the third plot of Fig. [2] is obtained using 
a different averaging procedure, and is therefore not a direct 
fit to the data shown. 

most terms in the aforementioned eigenvalue expansion 
will effectively remain constant. To account for the few 
modes which change during the decay, we tentatively 
write Mtrm (n>i) — c + c exp(a n/) + ci exp(ain/) + 
c 2 exp(a 2 n/) + . . . where dj are (real and negative) eigen- 
values, all of order one or smaller. Without loss of gener- 
ality, we can assume that for realistic values of m, a^ni 
remains sufficiently small to justify a linear expansion of 
the first exponential, leading to 

M TRM (ni) = c'+c aon/+ci exp(ain/)+c 2 exp(a 2 n/)+. . . 

(4) 

Averaging over nj yields 

M TRM (t,t w ) = c'+c aoaln(-^)+ Cl (^-) Al +c 2 (-^) A2 + . . . 

(5) 

where Aj = — a(l — e ai ) < and i — 1,2. Differentiating 
with respect to t, and re-arranging the terms, we finally 
obtain 

t r TRM {t) - c Q a a = Cl \i(—) Xl + c 2 A 2 (-^) A2 . (6) 

The above expression has the same functional form as 
Eq. [I] with a = coaoa, b\ — c\\\ and 6 2 = c 2 A 2 . The 
weak T dependence of both exponents and pre-factors 
should be expected as quakes are exothermal [19]. The 



near T independence of a implies, banning unlikely can- 
cellations, that a is itself temperature independent, pre- 
cisely as required by record dynamics. 

All physical observables are simply related in the lin- 
ear response regime [IfJ. E.g. the 'relaxation rate', 
which in the present notation is given by S(t) — 
— 8Mtrm / d\og(t — tw), can be cast into the form 
S(t obs ) = -t obs r TRM (t bs + t w ), where, as we recall, 
t Q ft s = t — t w . Inserting Eq. [5] for Ttrm an d plotting 
the outcome versus log 10 t b s reproduces the well known 
shape of S(t): a broad maximum is present at t f, s = t w , 
and a flat asymptotic value, equal to a, is reached for 
t bs ^ t w . The asymptotically logarithmic TRM de- 
cay and, equivalently, the constant asymptotic value of 
S are widely observed in complex systems, To name 
a few cases, they are seen in switchable mirrors after 
UV illumination [27l|. in the field-cooled magnetization 
of spin glasse s I 28l| . and in the magnetic creep of the 
ROM model 12911 of magnetic flux creep in type 11 su- 
perconductors [211 ]. In the asymptotic regime, the rate 
is Txrm oc 1/t, and the time at t w at which the pertur- 
bation is switched off is thus forgotten. Importantly, the 
energy decay rate of the EA spin glass [l5| and of a p-spin 
model [lij is also (nearly) proportional the reciprocal of 
the age, r^(t) oc 1/t. (Recall that, inconveniently, the 
symbol t w is used in Refs. [3, [3] in lieu oft.) In conclu- 
sion, relaxation in response to the initial thermal quench, 
which is never forgotten, seems to emerge as the main 
physical process during non-equilibrium aging, whether 
or not an added field is present. This is fully consistent 
with the subordination assumption used to interpret the 
present data. 

SUMMARY AND CONCLUSIONS 

Using available spin-glass thermoremanent magnetiza- 
tion data and a general method of analysis, we have ex- 
tracted the intermittent properties of the TRM decay and 
interpreted them theoretically using record dynamics. In 
combination with a previous investigation of the configu- 
ration autocorrelation function of the Edwards- Anderson 
spin-glass p^j . the present results suggests that conju- 
gate autocorrelation and response functions may inherit 
significant statistical properties from the quakes. This 
can lead to Fluctuation Dissipation-like relations [13] be- 
tween the two, even in an out-of-equilibrium situation, 
once the transient effects of the field switch have died 
out. 

The record dynamics approach melds real space fea- 
tures, i.e. the independent thermalized domains where 
quakes are initiated, with configuration space features, 
i.e. the scale invariance of the energy landscape associated 
to each domain, a concept repeatedly stressed in hierar- 
chical models [g, 0, S] ■ Similar ideas have been applied 
to 'non-thermal' dynamics, e.g. memory effects in driven 
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dissipative systems are linked to the marginal stability 
of the attractors selected by the dynamics [H, [U |32| 
and can also be understood by record dynamics. By con- 
struction, the description becomes invalid near the final 
equilibration time, since the postulate irreversibility of 
quakes cannot be maintained. 

The complex interplay between external noise and drift 
in glassy dynamics is a pivotal issue in non-equilibrium 
statistical physics, and many of its aspects are still only 
partially understood. Further insights could be obtained 
by applying the present method to linear response func- 
tions in other situations, e.g. aging with a small temper- 
ature step. The intermittency of the heat loss has been 
studied in some detail in this situation [5(, confirming 
that, in the Edwards- Anderson model, the largest energy 
barrier overcome in the past evolution sets the time scale 
for future intermittent events. One recent experimental 
finding is that the TRM signal loses its t w dependence 
in the extreme limit t w ^> t [33|]. As proposed in that 
paper, the 'post-aging' decay is intrinsically related to 
the same mechanisms producing the usual aging effects, 
and is related to the memory of the initial state set up 
by the cooling procedure. We expect that experimental 
investigations of intermittency in mesoscopic-scaled sys- 
tems will further elucidate this and other memory effects 
in complex dynamics. 
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Caption for Fig. 1 

(Color on line) (a): The circles show the PDF of the 
magnetization change SM for isothermal aging at T = 
0.83T g , with the field cut at t w = 100s. The values of 
SM are taken over short intervals St = 1.045s within 
the observation interval [1000, 4000]. Diamonds show the 
PDF of SM from the same data stream, but now taken 
over larger intervals St = 3 x 1.045s. This leads to a 
much stronger intermittent component on the left wing 
of the PDF. The full lines are least square fits of the 
positive values of SM to a Gaussian with zero average. 
(b): From the same data, the average and the variance 
(the latter scaled as shown) of SM 's are plotted versus 
St. The dotted line is obtained by a least square error fit. 

Caption for Fig. 2 

(Color on line) Panels 1-5: the main plots show the 
estimated values of t Ttrm — a versus t/t w , for t w = 50 
(right pointing triangles, unavailable for T = 0.6T g and 
T = 0.9T g ), 100 (circles), 300 (squares), 630 (diamonds), 
1000 (pentagrams), 3600 (hexagrams), 6310 (asterisks) 
and (10000) (left pointing triangles). The lines are given 
by Eq. [TJ The inserts compare the TRM decay measured 
at t w — 100s (red line), with the theoretical estimates 
(blue circles) obtained by integrating the fitted decay 
rate, see Eq. [2] The abscissa t/t w — 1 is the ratio of 
the observation time to the field removal age. In panel 
6, the parameter a and the two exponents Ai and A2 are 
plotted vs. temperature. The lines are guides to the eye. 

Caption for Fig. 3 

(Color on line) The main plot shows the unnormalized 
PDF (circles) of the quantity SM/St — a. The values of 
SM are collected from data streams taken with different 
t w and with T — 0.83T g . The denominator St is the re- 
peat time in each measurement, All data are taken within 
intervals [t, 1.254], where the value of t changes with t w , 
keeping a constant ratio t/t w = 1.13. The parabola and 
the straight line are fits to a Gaussian and an exponen- 
tial, respectively. The insert shows the t/t w dependence 
of t rTRAft — a (dots), where ttrm is estimated by aver- 
aging SM I St over the ensemble of data points available 
at each t/t w . The full line, which is lifted from the third 
plot of Fig. [21 is obtained using a different averaging pro- 
cedure, and is therefore not a direct fit to the data shown. 



